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Abstract 

The propagation of high-energy cosmic rays in the Galaxy is investigated. Solutions 
of a diffusion model are combined with numerically calculated trajectories of parti- 
cles. The resulting escape path length and interaction path length are presented and 
energy spectra obtained at Earth are discussed. It is shown that the energy spec- 
tra for heavy elements should be flatter as compared to light ones due to nuclear 
interactions during the propagation process. The obtained propagation properties 
of ultra-heavy elements indicate that these elements could play an important role 
for the explanation of the second knee in the cosmic-ray energy spectrum around 
400 PeV. 
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1 Introduction 



The explanation of the origin of super-high 
energy cosmic rays is one of the unsolved 
problems in astrophysics. The energy spec- 
tra at the sources are not identical to the 
observed spectra at Earth. Hence, studying 
the sources is closely related to investiga- 
tions of cosmic-ray propagation processes 
in the Galaxy. For the latter, a detailed 
knowledge of the structure of the magnetic 
fields is important. Unfortunately, the con- 



figuration of the galactic magnetic field re- 
mains an open question — different models 
can explain the experimental data [1,2,3,4]. 

How cosmic rays are accelerated to ex- 
tremely high energies is another unsolved 
problem. Although the popular model of 
cosmic-ray acceleration by shock waves in 
the expanding shells of supernovae (see 
e.g. [5,6,7]) is almost accepted as "stan- 
dard theory," there are still serious diffi- 
culties. Furthermore, the question about 
other acceleration mechanisms is not quite 
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solved, and such mechanisms could lead to 
different cosmic-ray energy spectra at the 
sources [1]. 



Different concepts can be verified, calcu- 
lating the primary cosmic-ray energy spec- 
trum, making assumptions on the density 
of cosmic-ray sources, the energy spectrum 
at the sources, and the configuration of 
the galactic magnetic fields. The diffusion 
model may be used in the energy range 
E < 10^'^ eV, where the energy spectrum 
is obtained using the diffusion equation for 
the density of cosmic rays in the Galaxy. 
At higher energies this model ceases to be 
valid, and it becomes necessary to carry 
out numerical calculations of particle tra- 
jectories for the propagation in the mag- 
netic fields. This method works best for 
the highest energy particles, since the time 
required for the calculations is inversely 
proportional to the particle energy. 



2 Assumptions 

High isotropy and a comparatively long re- 
tention of cosmic rays in the Galaxy (~ 10^ 
years for the disk model) reveal the diffu- 
sion nature of particle motion in the inter- 
stellar magnetic fields. This process is de- 
scribed by a corresponding diffusion tensor 
[1,3,8]. The steady-state diffusion equation 
for the cosmic-ray density N{r) is (neglect- 
ing nuclear interactions and energy losses) 

-WiD,j{r)WjN{r)^Q{r). (1) 

Q{r) is the cosmic-ray source term and 
Dij{r) the diffusion tensor. 

Under the assumption of azimuthal sym- 
metry and taking into account the predom- 
inance of the toroidal component of the 
magnetic field, eq. (1) is presented in cylin- 
drical coordinates as 



Therefore, a calculation of the cosmic-ray 
spectrum in the energy range 10^^ — 10^^ eV 
has been performed in a combined ap- 
proach: solutions of a diffusion model are 
used at low energies and particle trajec- 
tories are numerically integrated at high 
energies. 



In Sect. 2 the basic assumptions for the 
diffusion model will be described. The re- 
sults obtained with the propagation model 
are presented in the subsequent sections. 
The calculated propagation path length 
and interaction probability of cosmic rays 
will be discussed in Sect. 3 and Sect. 4, re- 
spectively. Finally, the energy spectra are 
presented in Sect. 5, followed by a discus- 
sion of the results (Sect. 6). 



1 d ^ d d ^ d d ^ d 
r dr dr dz dz dz dr 



r dr dz 



N{r,z)^Q{r,z){2) 



where N{r, z) is the cosmic-ray density 
averaged over the large-scale fluctuations 
with a characteristic scale L ~ 100 pc 
[3]. D± DC E"^ is the diffusion coefficient, 
where m is much less than one (m ~ 0.2), 
and Da oc E the Hall diffusion coefficient. 
The infiuence of Hall diffusion becomes 
predominant at high energies (> 10^^ eV). 
The sharp enhancement of the diffusion 
coefficient leads to an excessive cosmic- 
ray leakage from the Galaxy at energies 
E > 10^'^ eV. To investigate the cosmic-ray 
propagation at such energies it becomes 
necessary to calculate the trajectories for 
individual particles. 
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Such a numerical calculation of trajectories 
is based on the solution of the equation of 
motion for a charged particle in the mag- 
netic field. The calculation was carried out 
using a fourth-order Runge-Kutta method. 
Trajectories of cosmic rays were calculated 
until they left the Galaxy. Testing the dif- 
ferential scheme used, it was found that 
the accuracy of the obtained trajectories 
for protons with an energy E — 10^^ eV af- 
ter passing a distance of 1 pc amounts to 
5 • 10~^ pc. The retention time of a proton 
with such an energy averages to about 10 
million years, hence, the total error for the 
trajectory approximation by the differen- 
tial scheme is about 0.5 pc. 

The magnetic field of the Galaxy consists 
of a large-scale regular and a chaotic, irreg- 

— * — * — * 

ular component B — Br eg + B^rr- A purely 
azimuthal magnetic field was assumed for 
the regular field 



B^ = 0, B,f = 0, 
S<^ = l/xG-exp , (3) 

\ ^0 ^0 J 

where ZQ — h kpc and Tq = 10 kpc are con- 
stants [3]. These values are adopted from 
Ref. [3] to ensure the same conditions for 
both methods, i.e. trajectory calculations 
and the diffusion approach. The irregular 
field was constructed according to an algo- 
rithm used in [9], that takes into account 
the correlation of the magnetic field inten- 
sities in adjacent cells. The radius of the 
Galaxy is assumed to be 15 kpc and the 
galactic disk has a half-thickness of 200 pc. 
The position of the Solar system was de- 
fined at r = 8.5 kpc, (j) = 0°, and z = kpc. 
A radial distribution of supcrnovae rem- 
nants along the galactic disk was consid- 
ered as sources [10]. 
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Fig. 1. Path length in the Galaxy for protons. 

The values for the diffusion model (Xdif) 
indicated by the soHd line. They are extrapo- 
lated to lower energies by the dashed dotted 
line. Also shown are predictions of a leaky-box 
model (A;;,, eq. (4)), a residual path length 
model {Xrp, cq. (5)), and an upper limit for a 
residual path length model according to the 
TRACER experiment (Atr) [H]. The hori- 
zontal line indicates the matter to be passed 
along a straight line from the galactic center 
to the solar system (Age)- 

3 Propagation Path Length 

Assuming an interstellar matter density 
rid — 1 cm~^ for the galactic disk and 
Ufi — 0.01 cm~^ for the halo, following 
Chap. 3 in Ref. [1], trajectory calculations 
were performed at energies above 0.1 PeV. 
The dependence of the path length on en- 
ergy was obtained from the dependence 
of the transport time for protons in the 
galactic disk and the halo. The resulting 
escape path length for protons as func- 
tion of energy is presented in Fig. 1 as 
Xdif. For heavier nuclei with charge Z the 
path length scales with rigidity, i.e. is re- 
lated to the values for protons Xp{E) as 
X{E,Z) - XP{E/Z). 

For protons at 4 PeV, the amount of tra- 
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versed material is approximately 0. 7 g/ cm^ . 
At higher energies, the calculated path 
length decreases as oc Between 
0.1 and 1 PeV the calculations yield a 
behaviour A oc E'^ with 6 = 0.2. The 
dashed dotted line indicates a trend be- 
low 0.1 PcV extrapolating the calculated 
values to lower energies using the slope 
obtained. This yields a path length around 
15 g/cm^ at 1 GeV, well compatible with 
measured values, see below. The relatively 
small S = 0.2 is much lower than the value 
usually assumed in Leaky-Box models 
{S 0.6). On the other hand, such a slope 
oiS — 0.2 follows from the diffusion model 
proposed in Ref . [3] . 

Measurements of the ratio of secondary to 
primary cosmic-ray nuclei at GeV energies 
are successfully described using a leaky-box 
model. For example, assuming the escape 
path length for particles with rigidity R and 
velocity /3 = f/c as 



A 



lb 



26.7/3 g/cm^ 



\1.0 GYj 



^ \1A GV J 



-1.4 



(4) 



various secondary to primary ratios, like 
boron to carbon, phosphorus to sulfur, 
and sub-iron to iron, as obtained by the 
ACE/CRIS and HEAO-3 experiments, can 
be described consistently in the energy 
range from ~ 70 MeV to ~ 30 GeV [12]. 
A similar approach is the residual path 
length model [13], assuming the relation 



A, 



rp 



6.0 



R 



10 GV 



-0.6 



+ 0.013 



cm^ 



K5) 



for the escape path length. Recent mea- 
surements of the TRACER experiment 
yield an upper limit for the constant 
term in the residual path length model of 
0.15 g/cm^ [11]. The three examples are 



compared to the predictions of the diffu- 
sion model in Fig. 1. 

Extrapolating these relations to higher 
energies leads to extremely small values 
at PeV energies due to the strong de- 
pendence of the path length on energy 
(oc E~^'^). For the approaches according 
to eq. (4) and eq. (5) above 10^ GeV the 
traversed matter would be less than the 
matter passed along a straight line from 
the galactic center to the solar system 
Xgc = 8 kpc • 1 proton/cm^ 0.04 g/cm^. 
This value is indicated in Fig. 1 as dot- 
ted line. For the upper limit obtained by 
TRACER the path length at 10^ GeV 
amounts to about 0.17 g/cm^. This value 
does not exclude diffusion, as the larmor 
radius vl = 1.08 pc • E[PeV]/{Z ■ B[nG]) 
of a proton with an energy of 0.1 PeV in 
the galactic magnetic field (5 = 1 //G) is 
about 0.1 pc and 0.17 g/cm^ correspond 
to about 34 pc or SAOtl at 0.1 PeV. The 
value may be even larger when consider- 
ing cosmic-ray propagation in the galactic 
halo. However, in this approach the cosmic- 
ray sources are expected to be closer to the 
solar system as in the diffusion model, since 
the latter yields a path length of about 
1.5 g/cm^ at 0.1 PeV. This value follows 
from the assumption that the sources are 
disributed like supernova remnants and, 
hence, the question may arise whether this 
assumption is true. 

At low energies (GeV regime) the ratio of 
primary to secondary nuclei is measured 
to derive information about cosmic-ray 
propagation. At higher energies in the PeV 
domain such information is experimen- 
tally not accessible and other information 
is used to set constraints on propagation 
models. One possibility is the anisotropy 
of the arrival direction of cosmic rays. It 
can be characterized by Rayleigh ampli- 
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Fig. 2. Rayleigh amplitudes as function 
of energy for various experiments, taken 
from Ref. [14]. The results obtained by Mt. 
Norikura [15], EAS-TOP [16,17], Akeno [18], 
Adelaide [19], and KASCADE [14] are com- 
pared to predictions of leaky-box models [20] 
and a diffusion model [21]. For the diffusion 
model, predictions for primary protons, iron 
nuclei, and all particles are displayed. 
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Fig. 3. Interaction length as function of energy 
for different elements based on cross-sections 
according to the hadronic interaction model 
QGSJET. The nuclear charge numbers are in- 
dicated. 

the Galaxy at PeV energies. 



tudes. Values observed by different exper- 
iments are compiled in Fig. 2 [14]. They 
are compared to predictions of propaga- 
tion models. Two versions of a leaky-box 
model [20], with and without reaccelera- 
tion are represented in the figure. Leaky- 
box models, with their extremely steep 
decrease of the path length as function of 
energy (A oc E~^'^), yield relative large 
anisotropics even at modest energies be- 
low 10® GeV, which seem to be ruled out 
by the measurements, see also [20]. The 
measured values are almost an order of 
magnitude smaller. On the other hand, a 
diffusion model [21], which is based on the 
same basic idea [3] as the present work, 
predicts relatively small values at low en- 
ergies and a modest rise only. In the figure, 
predictions for pure protons and iron nu- 
clei are shown together with calculations 
for a mixed composition. The predicted 
Rayleigh amplitudes are compatible with 
the measured values. This may indicate, 
that diffusion models are a more realistic 
description of cosmic-ray propagation in 



4 Interaction Probability 

To calculate the interaction length of nuclei 
in the interstellar medium, nuclear cross 
sections according to the hadronic interac- 
tion model QGSJET [22] have been used. 
The cross sections for nuclei with mass 
number A > 1 have been parameterized 
using 

a{E)[mh] = a{E)A'^^^\ (6) 
where a and (3 are approximated as 

a{E) = 50.44 - 7.93 lg(^^^ 

+0.611g^(^) (7) 

and 

= 0.97- 0.022 Ig(^). (8) 
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Mass number In A 

Fig. 4. Parameters Aq at 10 GeV (filled sym- 
bols, left scale) and e (open symbols, right 
scale) to describe the interaction length as 
function of energy according to eq. (9) for nu- 
clei with mass number A. 

The corresponding interaction lengths for 
four different elements are presented in 
Fig. 3. The values decrease slightly as func- 
tion of energy. Values for protons are in 
the range 55 — 20 g/cm^, the values de- 
crease as function of nuclear charge and 
reach values < 1 g/cm^ at all energies for 
the heaviest elements. 

The dependence of the interaction length 
on energy has been approximated by the 
relation 



A 



int 



(9) 



The corresponding values for all elements 
are depicted in Fig. 4 as function of the 
logarithm of the mass number A. The val- 
ues for Ao at 10 GeV decrease from about 
50 g/cm^ for hydrogen to about 0.9 g/cw? 
for uranium. The small energy dependence 
is reflected in the values for e ranging from 
0.05 for protons to 0.01 for the heaviest 
element. The lines illustrate the approxi- 
mately hnear dependence of Ig Aq and e on 
In A. Small deviations from this behavior 
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Fig. 5. Fraction of nuclei surviving without in- 
teraction in the Galaxy for different elements 
as function of energy. 

are visible for the light elements helium, 
lithium, beryllium, and boron. The small 
values for e motivate why Ai„f is frequently 
assumed to be constant. 

Using the path length Xdif and the interac- 
tion length obtained, the interaction prob- 
ability for different nuclei has been calcu- 
lated. Nuclear fragmentation is taken into 
account in an approximate approach [23]. 
It should be pointed out that a nuclear frag- 
ment conserves the trajectory direction of 
its parent if Z/A in question is the same 
as for the primary nucleus and for most 
stable nuclei the ratio Z/A is close to 1/2. 
The resulting fraction of nuclei which sur- 
vive without an interaction is presented in 
Fig. 5 for selected elements. 

In the energy range shown protons reach 
the Earth almost undisturbed without any 
interaction, less than 5% of them interact. 
On the other hand, for heavier elements 
nuclear interactions play an important role 
and a significant fraction of them interacts 
at low energies. However, it turns out that 
even for ultra-heavy elements (heavier than 
iron) at the expected knees (at an energy 
of approximately Z ■ 4.5 PeV), see below. 
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Fig. 6. Calculated spectra of protons for the 
diffusion model (solid line) and the numerical 
trajectory calculations (dotted line). 

more than about 50% of the nuclei survive 
without interactions. This is an important 
result, which could help to understand the 
origin of the second knee in the all-particle 
spectrum around 400 PeV, as shall be dis- 
cussed below. It should be noted that the 
fraction of surviving nuclei is even larger 
for a leaky-box model, with its low path 
length at such high energies. 



5 Energy Spectra 



The results for the calculations of the 
cosmic-ray proton spectrum are presented 
in Fig. 6. They were obtained using the 
diffusion model and numerical calculations 
of trajectories. It is evident from the graph 
that both methods give identical results 
up to about 3 ■ 10^ GeV. At higher ener- 
gies there is a continuous decrease of the 
intensity in the diffusion spectrum, which 
has its origin in the excessive increase in 
the diffusion coefficient that leads to large 
leakage of particles from the Galaxy. 

An energy of 10® GeV can be taken as the 
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Fig. 7. Proton flux as obtained from various 
measurements, for references see [28], com- 
pared to the spectra shown in Fig. 6 (black 
lines) and the poly-gonato model [26] (grey, 
dashed line). 

conventional boundary to apply the diffu- 
sion model. At this energy the results ob- 
tained with the two methods differ by a fac- 
tor of two and for higher energies the dif- 
fusion approximation becomes invalid. 

Although the knee in the all-particle spec- 
trum has been observed more than 40 
years ago [24] , it was only recently that ex- 
perimental spectra for groups of elements 
became available. The KASCADE air 
shower experiment derived energy spectra 
for five groups of elements, namely pro- 
tons, helium, CNO, silicon group, and iron 
group [25]. The spectra exhibit a fall-off 
for individual elements at high energies. 
These results and the data available from 
other experiments are compatible with the 
poly-gonato model [26], assuming a knee 
for each element at an energy of about 
Z -4.5 PeV [27]. 

In the following, we compare the predicted 
spectra already shown in Fig. 6 to direct 
and indirect measurements of the primary 
proton flux in Fig. 7. The predicted spec- 
tra are normalized to average experimental 
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values at 1 PeV. In the range depicted, al- 
most no difference is seen between the two 
approaches. The relatively steep decrease 
of the measured flux at energies exceeding 
4 PeV is not reflected by the predictions. 
On the other hand, the data are described 
reasonably well by the poly-gonato model 
[26], shown in the figure as well. The ob- 
served change in the spectral index A7 f« 
2.1 according to the poly-gonato model has 
to be compared to the value predicted by 
the diffusion model. In the latter the change 
should be 1 — m ~ 0.8 [3]. The observed 
value is obviously larger, which implies that 
the remaining change of the spectral shape 
should be caused by a change of the spec- 
trum at the source, e.g. due to the max- 
imum energy attained in the acceleration 
process. 

The maximum energy and, therefore, the 
energy at which the spectrum steepens 
depends on the intensity of the magnetic 
fields in the acceleration zone and on a 
number of assumptions for the feedback of 
cosmic rays to the shock front. The uncer- 
tainty of the parameters yields variations 
in the maximum energy predicted by dif- 
ferent models up to a factor of 100 [6,29]. 
It seems, there is no consensus about what 
the "standard model" is considered to pre- 
dict. For the time being, it is difficult to 
draw definite conclusions from the compar- 
ison between the experimental spectra for 
different elemental groups and the "stan- 
dard model" of cosmic-ray acceleration at 
ultra high energies. 



6 Discussion 

The energy spectra for individual elements 
measured at the Earth with GeV and TeV 
energies can be described by power laws 



dN/dE oc E'^ with values for the spectral 
index 7 in the range -2.46 to -2.95 for ele- 
ments from hydrogen to nickel [30,26]. The 
measurements seem to indicate that the 
steepness of the energy spectra at Earth 
depends on the mass of the nuclei, heavier 
elements seem to have flatter spectra. At 
higher energies in the PeV domain the mea- 
sured spectra are compatible with the as- 
sumption of a knee for individual elements 
at about 4.5 PeV [26,27]. 

The energy spectrum of cosmic rays at their 
source Q{E) is related to the observed val- 
ues at Earth N{E) as 

with the escape path length X^sc and the in- 
teraction length Xint- Values for the former 
are presented in Fig. 1 and for the latter 
in Fig. 3. The relation between N{E) and 
Q{E) is governed by the absolute values 
of Xesc and Xint as well as their respective 
energy dependence. The interaction length 
Xint is almost independent of energy, the 
values for e in eq. (9) are smaller than 0.05, 
see Fig. 4. On the other hand, the propaga- 
tion path length Xesc decreases as function 
of energy as Xesc E~^, with values be- 
tween S = 0.6 for Leaky Box models and 
S = 0.2 for the diffusion model described 
in this work (see Sect. 3). 

In the " standard picture" of galactic cosmic 
rays usually Xint > Xesc is assumed with an 
energy independent interaction length and 
an escape path length Xesc = oc E~^-^. 
Thus, the relation between N{E) and Q{E) 
is dominated by the energy dependence of 
the escape path length. To explain the spec- 
trum observed at Earth N{E) oc E^'^''^ the 
spectrum at the sources has to be Q{E) oc 
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Recent measurements by the HESS exper- 
iment of the TeV gamma ray flux from the 
shell type supernova remnant RX J1713.7- 
3946 yield a spectral index 7 = — 2.19 ± 
0.09 ± 0.15 [31]. For the first time, spec- 
tra have been derived directly at a poten- 
tial cosmic-ray source. The results are com- 
patible with a nonlinear kinetic theory of 
cosmic-ray acceleration in supernova rem- 
nants and imply that this supernova rem- 
nant is an effective source of nuclear cosmic 
rays, where about 10% of the mechanical 
explosion energy are converted into nuclear 
cosmic rays [32] . 

These investigations show that the spec- 
trum at the source is Q{E) cx E'"^-^, a value 
relatively close to the "standard model". 
The escape path length required to obtain 
the spectra observed at Earth should be 
(X E^^'^, i.e. S is 0.1 smaller than in the 
" standard picture" . Two questions remain 
open: Do all potential cosmic-ray sources 
exhibit a similar spectrum? And, what frac- 
tion of galactic cosmic rays is accelerated in 
supernova remnants with spectra like RX 
J1713? 

Considering the spectra of nuclei observed 
at Earth, the second term in eq. (10) has to 
be taken into account. The vahics of Xint for 
protons are at all energies larger than the 
escape path length Xdif- Hence, for protons 
the escape from the Galaxy is the domi- 
nant process influencing the shape of the 
observed energy spectrum. For iron nuclei 
at low energies hadronic interactions are 
dominating {Xint ~ 2.5 g/cm^ < Xdif) and 
leakage from the Galaxy becomes impor- 
tant only at energies approaching the iron 
knee (26 • Ep). Finally, for elements heav- 
ier than iron, the interaction path length is 
smaller than the escape path length at all 
energies, except above the respective knees. 
For these elements the propagation path 



length exceeds the interaction path length 
by about an order of magnitude at low en- 
ergies. Hence, nuclear interaction processes 
are dominant for the shape of the observed 
spectrum. Taking this effect into account, 
it is expected that the energy spectra for 
heavy elements should be flatter as com- 
pared to light nuclei. The spectral index 
for iron nuclei should be about 0.13 smaller 
than the corresponding value for protons. 
Direct measurements indicate indeed that 
the spectra of light elements are flatter as 
compared to heavy elements [26], e.g. the 
values for protons 7p = 2.71 and iron ^ype — 
2.59 differ as expected. This more or less 
obvious effect has been pointed out earlier 
in several articles, see e.g. [33,34]. However, 
in the present work we extend the consid- 
erations to a wider energy range and ultra- 
heavy nuclei using contemporary estimates 
of the cross sections. 

The poly-gonato model extrapolates re- 
sults from direct measurements in the GeV 
and TeV regime to higher energies [26] . In 
the model fall-offs for individual elements 
proportional to the nuclear charge are as- 
sumed at energies of Ez — Z ■ Ep, with 
the value Ep — 4.5 PeV for protons. The 
knee is caused by the subsequential fall- 
offs for all elements, starting with protons. 
The second knee in the all-particle spec- 
trum around 400 PeV is interpreted as the 
end of the galactic component, where the 
heaviest elements (actinides) reach their 
fall-off energies (92 • Ep ^ 414 PeV). At 
this energy, elements heavier than iron are 
expected to contribute with about 40% to 
the all-particle flux [35]. The model uses an 
empirical relation for the spectral index. 
The spectral slope decreases as function of 
nuclear charge. In the following, this be- 
havior is investigated quantitatively with 
the propagation model discussed above. 
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Fig. 8. Energy spectra at Earth for elements with nuclear charge Z as indicated. The dashed 
lines represent spectra according to the poly-gonato model, the solid lines are expected from 
the diffusion model discussed, see text. Two solid lines are shown in each panel, representing an 
estimate for the uncertainties. 



As source composition, the abundances 
of elements from hydrogen to uranium as 
measured in the solar system [36] have been 
weighted with Z^'^. This choice is arbi- 
trary to a certain extent, but may be moti- 
vated by a higher efficiency in the injection 
or acceleration processes for nuclei with 
high charge numbers. The abundances are 
scaled with a factor which is identical for 
all elements to obtain approximately the 
absolute values as expected at the Earth 
according to the poly-gonato model. At the 
source, a power law oc has been as- 

sumed for all elements with a knee, caused 
e.g. by the maximum energy attained dur- 
ing the acceleration, at Z ■ 4.5 PeV, with 
a power law index —3.5 above the respec- 
tive knee. Using the derived propagation 
path length and interaction length, the 



amount of interacting particles has been 
determined. Secondary products gener- 
ated in spallation processes are taken into 
account, assuming that the energy per nu- 
cleon is conserved in these reactions. They 
are added to the corresponding spectra 
with smaller Z. The spectra thus obtained 
are compared to spectra according to the 
poly-gonato model in Fig. 8. 

Two features should be noted: The abso- 
lute fluxes at Earth are predicted quite 
well, especially when considering that only 
a simple scaling law has been introduced 
for the abundances at the sources, starting 
with the composition in the solar system. 
More important for the present discussion 
is the shape of the spectra. As expected, 
the shape of the proton spectrum is not 
influenced by the (few) interactions dur- 
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ing propagation and the difference of the 
spectral index at the source and at Earth 
7 = —2.71 [26] can be explained by the en- 
ergy dependence of the escape path length 
oc E~^-^. On the other hand, it can be rec- 
ognized that due to nuclear interactions 
the spectra for heavier elements are indeed 
flatter. The slopes obtained with the sim- 
ple approach for the CNO, silicon, and iron 
groups agree well with the steepness as 
expected from the poly-gonato model. For 
heavy elements at low energies secondary 
products generated in spallation processes 
play an important role for the shape of the 
spectrum. At low energies many nuclei in- 
teract due to the large escape path length 
and the small interaction length, thus, the 
spectra of nuclei without any interaction 
deviate from power laws. However, the 
spallation products of heavier elements at 
higher energies compensate the effect and 
the resulting spectra are again approxi- 
mately power laws, as can be seen in Fig. 8. 



7 Summciry 

The results obtained show the effectiveness 
of the combined method to calculate the 

cosmic-ray spectrum using a numerical cal- 
culation of trajectories and a diffusion ap- 
proximation. The propagation path length 
in the Galaxy and the energy spectra at 
Earth have been calculated to investigate 
the effect of escape from the Galaxy as pos- 
sible origin of the knee in the cosmic-ray 
energy spectrum. The calculated depen- 
dence of the propagation path length on 
energy suggests difficulties with the " stan- 
dard model" of galactic cosmic rays. An 
energy dependence of the propagation path 
length at low energies oc E~^-^ is expected, 
which requires a spectrum at the sources 



oc E^"^'^ in order to explain the observed 
spectra at Earth. To explain the relatively 
steep fall-off of the observed energy spec- 
tra for elemental groups at their respective 
knees, the modulation of the spectrum 
due to propagation solely is not sufficient. 
An additional steepening of the spectra at 
the source is necessary, e.g. caused by the 
maximum energy attained during acceler- 
ation. It can be concluded that the knee 
in the energy spectrum of cosmic rays has 
its origin most likely in both, acceleration 
and propagation processes. 

In the poly-gonato model the knee in the 
energy spectrum at Ep = 4.5 PeV is caused 
by a cut-off of the light elements and the 
spectrum above the knee is determined 
by the subsequent cut-offs of all heavier 
elements at energies proportional to their 
nuclear charge number. The second knee 
around 400 PeV ^ 92 ■ Ep is due to the 
cut-off of the heaviest elements in galactic 
cosmic rays. Considering the calculated 
escape path length and nuclear interac- 
tion length within the diffusion model, it 
seems to be reasonable that the spectra for 
heavy elements are flatter as compared to 
hght elements. The calculations show also 
that even for the heaviest elements at the 
respective knee energies more than about 
50% of the nuclei survive the propagation 
process without interactions. This may 
explain why ultra-heavy elements could 
contribute significantly (~ 40%) to the all- 
particle flux at energies around 400 PeV 
and thus explain the second knee in the 
energy spectrum. 
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